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We study radial wave packets produced by short-pulsed laser fields 
acting on Rydberg atoms, using analytical tools from supersymmetry- 
based quantum-defect theory We begin with a time-dependent per- 
turbative calculation for alkali-metal atoms, incorporating the atomic- 
excitation process. This provides insight into the general wave packet 
behavior and demonstrates agreement with conventional theory. We 
then obtain an alternative analytical description of a radial wave packet 
as a member of a particular family of squeezed states, which we call 
radial squeezed states. By construction, these have close to minimum 
uncertainty in the radial coordinates during the first pass through the 
outer apsidal point. The properties of radial squeezed states are inves- 
tigated, and they are shown to provide a description of certain aspects 
of Rydberg atoms excited by short-pulsed laser fields. We derive ex- 
pressions for the time evolution and the autocorrelation of the radial 
squeezed states, and we study numerically and analytically their be- 
havior in several alkali-metal atoms. Full and fractional revivals are 
observed. Comparisons show agreement with other theoretical results 
and with experiment. 
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I. INTRODUCTION 



When an atom is excited from its ground state by a short-pulsed laser field, a 
localized radial wave packet is produced that has behavior mimicking the classical 
radial motion of a charged particle in a Coulomb field 0, |[]. The properties and 
time evolution of such wave packets provide interesting experimental and theoretical 
opportunities to probe the interface between classical and quantum mechanics. 

Experiments have detected the periodic motion of the electron in a Rydberg atom 
excited by a short laser pulse, with period equal to the classical period of a particle in a 
keplerian orbit ^ |J. The motion of the wave packet is only partly classical, however. 
For excitations by a single laser pulse from the atomic ground state, the angular 
distribution is that of a p state. Also, the radial wave packet disperses over the course 
of time. However, after many Kepler times, the wave packet remnants eventually 
recombine into a packet that is close to the original shape and that oscillates radially 
with the keplerian orbital period. This recombined wave packet is called a full revival. 
In the intervals between the initial and the fully revived classical motions, subsidiary 
wave packets form. These are called fractional revivals, and they have orbital periods 
equal to rational fractions of the classical keplerian period. Both full and fractional 
revivals have been seen experimentally f7|, || [|, and descriptions of their behavior 
have been developed 

Several theoretical approaches have been used to study the properties of radial 
wave packets formed by excitation of Rydberg atoms with short laser pulses. Refs. 
|1|, § solve the Schrodinger equation numerically, while for gaussian shaped laser 
pulses and weak-field excitations, Refs. Q use perturbation theory. For a particular 
class of pulse shapes, Ref. is able to obtain nonperturbative solutions. Also, 
Ref. provides a theoretical description for the generation of Ramsey fringes with 
Rydberg wave packets. 

Since the radial wave packets are initially localized in the radial coordinates and 
partly follow the classical motion, a description in terms of some kind of coherent 
state O] might seem appropriate. However, standard coherent-state approaches to 



the hydrogenic Coulomb problem |16, 17, IS, 19, 20] and to short-pulse laser excitation 
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of Rydberg atoms |21. 22. 23] have yielded descriptions of the motion in a Coulomb 
potential with the electron moving on circular or elliptical orbits. None of these 
match the behavior of Rydberg atoms in short-pulsed laser fields with p-state angular 
distributions. Within this framework, the issue of an analytical construction of radial 
wave packets for Rydberg atoms prepared by excitation with short laser pulses with 
no external fields present has been an open problem. 

Recently, we have provided a framework for an analytical study of Rydberg wave 
packets |24j] and have discussed some of their properties for hydrogen. In the present 
work, we generalize this method to non-hydrogenic atoms, with particular emphasis 
on the alkali-metal atoms for which experiments have been performed. We also extend 
our analysis and obtain further results for the hydrogenic case. 

To incorporate non-hydrogenic features in the treatment, we address matters in 



the context of supersymmetry-based quantum-defect theory (SQDT) |25], ^6|. This 
theory describes the behavior of an excited Rydberg electron as that of a single 
particle in an effective central potential. The effective potential can be found by 
acting on the hydrogenic Coulomb potential with a symmetry transformation and 
adding a symmetry-breaking term. The symmetry transformation used is a quantum- 
mechanical supersymmetry [27]. The supersymmetry-breaking term incorporates 
nonhydrogenic contributions to the effective potential. Although the model is an- 
alytical, its combination of supersymmetry ideas with the notion of quantum defects 
P%| , |2"D| . pp| . |3~T| is such that exact energy eigenvalues are generated. The corresponding 
analytical eigenf unctions can be used to provide predictions for transition probabil- 
ities of alkali-metal atoms [26] and alkaline-earth-metal ions p2 , in good agreement 
with accepted values. The model also can be used to generate Stark maps for the 
alkali- metal atoms |33| that agree with experiment and standard theory |34|, [35|. The 
validity of the model breaks down inside the electronic core, although some bulk fea- 
tures of the fine structure of alkali-metal atoms are reproduced ||36|| . Extensions of 



these ideas apply to other situations including, for example, the Penning trap [p7 
For recent overviews of atomic supersymmetry, see Ref. ][8| . 

The analytical inclusion of quantum defects for Rydberg atoms suggests that 
SQDT could be a useful tool for the study of Rydberg atoms in short-pulsed laser 
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fields. We have identified two approaches as being especially promising in this regard. 
One is the inclusion of SQDT in conventional time-dependent perturbation theory to 
explore the formation and behavior of radial wave packets, in particular to show the 
presence of oscillatory motion between the orbital apsidal points. The results can 
then be compared with those obtained using numerical approaches [[1], 0] and WKB 
solutions in perturbation theory || 

The second approach we adopt is the use of SQDT in the construction of a family 
of analytical squeezed states, which we call radial squeezed states (RSS), that can 
represent radial wave packets formed in a Rydberg atom excited by a short laser pulse. 
For our purposes, given a hamiltonian and a commutation relation, a squeezed state 
can be taken as a wave function that satisfies the corresponding minimum-uncertainty 
relation at a particular time but with uncertainties not necessarily those of the ground 
state. For the simple harmonic oscillator, for example, the uncertainty product of a 
squeezed state has a sinusoidal dependence on time, with the uncertainties in coordi- 
nate and momentum space oscillating while the center of the wave packet follows the 
classical motion j39] . 

The analytical method we use to obtain the RSS is discussed in Sec. IV. The result 
is a three-parameter family of wave packets with the angular dependence of a p state. 
For an individual packet, the values of the three parameters are fixed in terms of ex- 
pectations of the energy, momentum, and position when the wave packet first reaches 
the outer apsidal point. This initialization condition agrees both with the results of 
other theoretical calculations and with the intuitive notion that the electron wave 
packet attains closest-to-classical behavior when localized at large quantum numbers 
immediately after its formation. The ensuing behavior is completely determined by 
the time-dependent Schrodinger equation. 

Since the RSS are given by a relatively simple analytical expression, their time 
evolution can be studied either analytically or numerically. We have used both ap- 
proaches. Quantities of interest are the time dependence of the RSS and of autocorre- 
lation functions with the initial packet. In what follows, we examine these for several 
different alkali-metal atoms. We treat issues such as the variation with the quantum 
defect of the orbital period and of the full or fractional revivals. 
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This paper is organized as follows. In the next section, we summarize some nec- 
essary background information on the salient features of SQDT and on the various 
approaches to the study of radial wave packets formed with short-pulsed laser fields 
applied to Rydberg atoms. In Sec. Ill, we use SQDT to perform a time-dependent 
perturbative calculation for the formation and time evolution of the radial wave pack- 
ets. The derivation and basic properties of the RSS are described in Sec. IV. Section 
V discusses their time-dependent properties, including their time evolution and their 
autocorrelation functions, comparing with other available theoretical results and with 
some experiments. A summary and our conclusions are given in Sec. VI. Throughout 
this work, we use atomic units with H — e — m e — 1. 

II. BACKGROUND 

In this section, we provide background material on atomic supersymmetry and on 
prior theoretical treatments of Rydberg atoms in short-pulsed laser fields. 

A. Atomic Supersymmetry 

Atomic supersymmetry involves a particular realization 7Z of the superalgebra 
sqm(2). The three generators of sqm(2) can be taken as the hamiltonian H ss of the 
super symmetric quantum system, the supersymmetry charge Q, and its conjugate 
Q^. They obey the relations 

[H ss ,Q] = [H ss ,Q^] = , {Q,Q^} = H SS . (1) 

In the realization 7Z, the hamiltonian for sqm(2) splits as the direct sum H ss = 
H + © H-. The component hamiltonians H + and H are one- variable differential 
operators satisfying the eigenequations 

d 2 

H ± y ±n = [-— + V±(x)}y± n = 6 n y± n , (2) 

where the potentials V + and V- are supersymmetric partners conventionally expressed 
in terms of ^-derivatives of a function U(x): 

V ± {x) = {\U i Yt\U" . (3) 
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The ground-state eigenvalue is zero and is associated with H + only. The spectra of 
the two hamiltonians are otherwise identical [fZ7] . 

When the hamiltonian H + is identified with the differential operator for the radial 
part of the hydrogen-atom Schrodinger equation expressed in spherical polar coordi- 



nates, the supersymmetry partner if_ can be constructed |25[ by fixing / and solving 
for V-i in Eq. (^). It can then be shown that the supersymmetric eigenfunction 
partners are R n i with n > 1 and R n ,i+i with n > 2. A physical interpretation of 
these results has been proposed |25j , leading to sets of nested supersymmetries inter- 
connecting eigenvalues and eigenstates of different atoms and ions. For example, in 
the exact symmetry limit where electron-electron interactions can be neglected, the 
s levels of lithium may be regarded as the supersymmetric partners of the hydrogen 
atom s levels. The s orbitals of sodium can then in turn be viewed as supersymmetric 
partners of the lithium s orbitals. Connections also exist among p and higher orbitals. 
An integer shift in the angular quantum number / accompanies each supersymmetry 
operation. 

In real atoms such as alkali metals, the exact supersymmetry is broken by cou- 
pling between the valence electron and the electronic core. One result is a shift in 
eigenenergies P5| , E n — > E n * = —l/2n* 2 , relative to hydrogenic values. The quantity 
n* = n — 5(n, I) incorporates these shifts via quantum defects 5(n, I). As n increases, 
the exact quantum defects rapidly approach asymptotic values 5(1). 

The basic idea of SQDT is to incorporate these shifts in an analytical one-particle 
hamiltonian ||26| . The model is obtained by adding to the exact supersymmetric 
hamiltonian specific supersymmetry-breaking terms, chosen so that the quantum- 
defect eigenenergies are reproduced while leaving the eigenfunctions analytical. A 
modified angular quantum number I* = I — 5(1) + 7(7) is introduced, where 1(1) is 
an integer playing the role of the supersymmetric shift. The SQDT hamiltonian is 
found via the replacement n, /, E n — > n*, /*, E n * in the radial Coulomb equation. The 
corresponding SQDT eigenfunctions are given by 



Tin* - I* 



T(n* + l* + 1) 



R^(r) = — 2 ^\ ; ^ ^ '(-) exp(--)L^ 1 (-) • W 



\ n *j \ n * / \n*' 



The full three-dimensional SQDT wave functions are therefore R n *i* (r)Yi m (8, 0). For 
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asymptotic quantum defects 5(1), these eigenfunctions form a complete and orthonor- 
mal set. Note that in the limit of vanishing quantum defects 5(1) and supersymmetry 
integers 1(1) the usual Coulomb eigenvalues and eigenfunctions are recovered. 

B. Radial Wave Packets 

The time-dependent Schrodinger equation for a hydrogen atom in a laser field is 

i^(t) = (H -fl-E)*(t) , (5) 

where Ho is the zero-field hamiltonian, \x is the atomic dipole moment, and E is the 
laser field. Consider an atom excited by a single photon from an initial state \i) to 
a range of Rydberg states by a laser pulse of length r and frequency lu. The electric 
field can be written as 

E(t) = £(t) ee~ lujt + c.c. . (6) 
Here, e is the polarization vector and £ (t) is a gaussian envelope for the laser pulse, 

e(t)=ft „p(_«) , (7) 

which is taken to be centered at t — with full width at half maximum equal to r. 
If we expand ^(t) in terms of hydrogenic wave functions, 

$>(t) = J2am m (t)e- iEnt R n iY lrn , (8) 

nlm 

and use the rotating wave approximation, we obtain equations for the time-dependent 
parameters a n i m (t). These equations may be solved numerically [|T], or perturba- 
tively for weak fields || |]. Quantities of interest, such as the probability distribution 
r 2 \^(r, t)\ 2 of the radial part of the wave function, may then be plotted as a function 
of time. Refs. |L[ f|, ]|J] use hydrogenic dipole moments and energy spacings for sim- 
plicity, although numerical methods incorporating standard quantum-defect theory 
can be applied. Refs. || [7], || contain numerical calculations for alkali-metal atoms. 

The results obtained by these procedures show that a radial wave packet forms 
shortly after the initial laser excitation. The packet subsequently oscillates between 
the radial apsidal points. Refs. [|l|, 0, |3|, |J] consider an 8 to 10 psec laser pulse 
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that excites a range of p states of hydrogen centered about the value n = 85 of 
n. The resulting wave packet is calculated to oscillate with a period equal to the 
classical Kepler period T c \ = 27m 3 ~ 93.3 psec. At the inner turning point, near 
the origin, the distribution r 2 |\I/(r, t)\ 2 of the packet is broad and exhibits many 
oscillations. As the packet returns to the outer turning point, near 2n 2 = 14450 a.u., 
the distribution becomes narrower and increases in amplitude, while the number of 
oscillations decreases significantly. At the first pass through the outer turning point, 
the wave packet is calculated to have close to minimum uncertainty in the radial 
coordinates Q, 0] . 

III. TIME-DEPENDENT PERTURBATION THEORY 

In this section, we show how to perform a time- dependent perturbative calculation 
that incorporates SQDT into the analysis. The calculation describes the formation 
and time evolution of a Rydberg atom excited from its ground state by a short laser 
pulse. The field is taken to have the form in Eq. @ with envelope £(t) as in Eq. (0). 

We begin by expanding the wave function describing alkali- metal Rydberg 
states in terms of the complete and orthonormal set of SQDT states R n *i*(r)Yi m (9, (p). 
For single photon excitations from the ground state, only p-state angular wave func- 
tions are needed. However, the procedure that follows can be generalized to multi- 
photon processes that excite d or higher states. We neglect any contributions from the 
continuum, which are negligible for the cases of interest. The expansion (§) therefore 
becomes 

tf(t) =Y w (6,<j ) )Y, a n(t)e- lE "* t Rn*i*(r) , (9) 

n 

where now n* = n — 5(1) and I* = 1 — 5(1) + 1(1), with 5(1) the p-state quantum 
defect and 1(1) the p-state supersymmetry integer. We then use the rotating wave 
approximation and solve for a n (t) to lowest order in perturbation theory. For times 
t 3> 0, we get 

a n = i(ri*,l*\j2- e\i)S(5 n ) , (10) 

where \i) denotes the ground state and S(5 n ) is the Fourier transform of S(t) with 
5 n = (E n * — Ei* — to). Note that the eigenenergies E n * = —l/2n* 2 , including the 
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ground-state energy incorporate the quantum defect. 

The dipole matrix elements (n*,l*\jl ■ e\i) can be evaluated analytically for dif- 
ferent alkali- metal atoms Equation (^) specifies the form of the radial wave 
functions R n *i*(r). Then, for fixed 9 and up to a normalization constant N, the 
radial probability density is given by 

r 2 |*(t)| 2 = iVr 2 | (n*,r|r|l,0> exp (-^^f^) e"^-* f . 

(11) 

In this expression, we have taken the laser frequency u as tuned to excite a range of 
energies centered on E h * , 

00 = E n * - Ei* . (12) 

Equation flTT| ) holds for t ^> 0. We can modify the procedure and evaluate r 2 |\l/(i) | 2 
for t = as well. The resulting calculation shows that the answer can be written by 
making the substitution 



-iE n *t 1 
C 2 



1 + erf W 161^2 )J (13) 



in Eq. (|TT|), where erf(z) is the error function. 

As a first example, consider hydrogen. This element has 5(1) = 1(1) = 0, so 
n* = n. For purposes of comparison with Ref. [0], we take n* = n = 85 and r = 8 
psec. Figure 1 displays the radial probability distribution f(r) = r 2 \^/(t)\ 2 at several 
different times. From left to right, these are t = |T d ~ 10.3 psec, |T cl ~ 20.7 psec, 
|T cl ~ 31.1 psec, and |T c i ~ 41.5 psec, where T c \ = 27m 3 ~ 93.3 psec is the classical 
orbital period. No vertical scale is shown because the packet is unnormalized. The 
figure uses Eq. (|TTJ) evaluated numerically over the subset of states n = 75-95, which 
we determined to be sufficient for high accuracy. 

The graph reveals that during formation the radial wave packet moves away from 
the origin, increasing in amplitude as it approaches the outer turning point r c \ w 2n 2 = 
14450 a.u. It then decreases in amplitude as it starts to bounce back towards the 
origin. Our Fig. 1 agrees with Figure 1 of Ref. || in its gross features, although there 
are differences in the detailed shapes of the packets. Ref. uses a WKB solution 
to evaluate the Coulomb functions, whereas we use Eq. (|TT|) directly and evaluate 
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the associated Laguerre polynomials in R n *i*(r) with high numerical accuracy for a 
subset of states. 

As a second example, we look at rubidium. This has 5(1) ~ 2.65 and 1(1) = 3. 
For purposes of comparison with the hydrogenic case, we again choose r = 8 psec and 
n = 85. This corresponds to n* ~ 82.35. As we show in Sec. IV, the inequivalence 
of n and n* means that the classical orbital period is reduced from the hydrogenic 
value to T c * = 27m* 3 ~ 84.9 psec. Figure 2 displays the radial probability distribution 
f(r) for rubidium at the same times as in Fig. 1. The overall behavior is similar to 
that of hydrogen. However, comparison of Figs. 1 and 2 shows that the outer turning 
point r*j for rubidium is slightly smaller than for hydrogen. This difference is due to 
the quantum defects. In fact, as we show in Sec. IV A, the outer turning point is 
r*j « 2n* 2 ~ 13560 a.u. here. 

IV. RADIAL ATOMIC SQUEEZED STATES 

The results in the previous section demonstrate that the analytical framework 
provided by SQDT is suitable for the perturbative study of the time evolution of 
radial wave packets for different alkali-metal atoms. However, the development and 
application of an analytical treatment of the problem that avoids perturbation theory 
is evidently of interest and is the principal goal of this paper. 

The uncertainty product ArAp r is close to its minimum value during the first 
passage of the wave packet through the outer turning point [p], |2|]. This suggests that 
we seek an analytical description of a radial wave packet at its outer turning point as 
some type of minimum-uncertainty wave function. 

In this section, we use SQDT to derive a family of analytical wave packets that 
obey minimum-uncertainty relations for certain special variables and that have un- 
certainty product ArAp r close to the minimum value. We call these packets radial 
squeezed states (RSS). The use of SQDT ensures that the RSS are relevant for all 
alkali-metal atoms, not merely hydrogen, and permits us to take into account the 
quantum defects of the different alkali-metal atoms used in experiments. Discussion 
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of the time evolution of the RSS is deferred to Sec. V. 



A. Classical Motion 

The construction of close-to-classical (approximate minimum-uncertainty) wave 
functions is aided by an understanding of classical behavior. In this subsection, we 
consider the motion of a classical particle in a central potential that can be viewed 
as the classical limit of the effective potential in SQDT. 

In general, the motion of a particle in a central potential U (r) is governed by the 
lagrangian 

L = \r 2 + \r 2 6 2 - U(r) . (14) 

The classical angular momentum / = r 2 9 is constant and so can be used to reduce 
the equations of motion to a single equation for the radial coordinate r. 

The SQDT uses an /-dependent effective central potential U(r). One perspective 
is therefore that it replaces the intractable many-body action for the multi-electron 
atom with an infinite set of single-particle effective actions, labeled by I. Each member 
of this set describes the behavior of the system for the given value of the angular 
momentum. In the classical case, I is a continuous variable, and the central potential 
U (r) can be taken as 

1 7*2 _ 72 

This potential generates an effective radial hamiltonian 



H *= l 2P 2 r + ^-l=E* , (16) 



where p r = r is the radial momentum and E* is the energy. Note that when I* = /, 
the classical SQDT reduces to the usual Coulomb case. 

For negative E*, the radial motion is oscillatory between two values r 1>2 of r 
corresponding to the apsidal points of the orbit and given by 



1- (l ± v/l-2|£*r) . (17) 
The classical orbital period T* x for the SQDT can be defined as the time taken to 



n ' 2= 2\E* 
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move from r\ to r 2 and back. Using Eq. (jig), a short calculation gives 

1 cl ~ 



2|£*| 3 

Bertrand's theorem implies that for the SQDT central potential (|l^) the orbits 
are not closed. However, the classical equations of motion can be solved exactly, as 
might be expected from the analytical nature of the quantum theory. The solution 
for the classical orbit can be found by eliminating t from Eq. (|T6|), using d6 = ldt/r 2 , 
and integrating. We thereby obtain the orbit equation 

i = i(l + e coB[f(0-0o)]) . (19) 

Here, 8q is a constant of the integration, and we have defined 

I* 

f = j (20) 

and 

e = y/l- 2\E*\l* 2 . (21) 

Equation ( |T9"| ) describes a precessing ellipse, with semimajor axis a = 1/2\E*\ and 
eccentricity e. The apsidal points of the orbit are r 12 = a(l ± e). The precession 
in one orbital period is A6 = 2tv(1 — 1/f)- Evidently, the precession is a classical 
consequence of the presence in the SQDT of the supersymmetry-type shifts and the 
quantum defects. Note that for different atomic levels the quantity / in the SQDT can 
be larger or smaller than one, so the precession may be clockwise or counterclockwise. 
Also, if I* = I then / = 1 and the precession disappears. 

For convenience in comparing with the quantum theory, it is useful to write E* as 

E* = . (22) 

2n* 2 v ' 

The choice of notation reflects that in the quantum theory n* is the quantized, shifted 
principal quantum number. However, at the classical level, it is a convenient, con- 
tinuous variable. The classical turning points and the classical orbital period can be 
expressed in terms of n* as 



r li2 = n* 2 l±Wl-_ (23) 
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and 

T c * = 2vm* 3 . (24) 

The orbital eccentricity becomes 

H 1 -^ • (25) 

For single-photon excitation of a Rydberg atom to a wave packet with large average 
principal quantum number, the ratio l*/n* is small. This means that the RSS we de- 
rive below and the radial wave packets experimentally constructed have corresponding 
classical orbits that are highly elliptical, with e ~ 1. 



B. Oscillator Description of Radial Motion 

We have found that direct attempts to construct minimum-uncertainty wave func- 
tions for the SQDT meet intractable difficulties. In fact, the same is true of the pure 



hydrogenic case, as has been known since the time of Schrodinger |fL6j . We sidestep 
the issue by replacing r and p r with a new set of classical variables, R and P, in terms 
of which the description of the motion is close to that of a harmonic oscillator. The 
new radial variable R is chosen to have a simple sinusoidal dependence on the angle 
9. The resulting equations are relatively easy to handle at the quantum level, so that 
minimum-uncertainty packets can be found analytically. This elegant technique was 
originally introduced in Ref. [H5] to construct 'minimum-uncertainty coherent states' 
for the Coulomb problem, which correspond in the present context to states with high 
angular momentum (and small eccentricity). 

For simplicity, we pick 9 such that f9 = |. The new radial variable is then given 

as 



a = ;-^f = • (26) 



Its time derivative is 



A=^cos/0 , (27) 



and the corresponding conjugate momentum P is 

r2 R ^ 



P = --R=-— cos f9 . (28) 
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The classical equations of motion can be written 



R = ~L P > ( 29 ) 



r 

The new classical variables satisfy the relations 



l 2 f 

P=-Lr . (30) 



r Z*rV 2Er 2 + 2r - I* 2 ' 1 ' 

p=1 f = 7 ' (32) 

These can be used to reexpress Eq. ( [16|) as 

^ 2 + ^ 2 = ^ • (33) 

This equation has the form of the energy equation for a simple harmonic oscillator of 
frequency I and energy e 2 /2f 2 . 

C. Derivation of Radial Squeezed States 

We next pass to the quantum level. The quantum radial hamiltonian for SQDT 



is 



H - U d2 + 2d ) 1 _|_ + X ) ( o 1; 

2 ^dr 2 rdr' r 2r 2 

The canonical radial coordinate r and its conjugate momentum p r = —i(d r + ^) obey 
the commutation relation [r,p r ] = i. The eigensolutions for this hamiltonian are given 
in Eq. (|). 

The quantum operators corresponding to the new classical variables R, P are 

R= r "/*(/* + 1) ' (35) 

P = P i = -- f {dr + -) ■ (36) 

f f r 



They obey the commutation relation 

i 1 



[ R i P \ = ---2 ■ (37) 
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The presence of the factor 1/r 2 in this equation reflects the unconventional choice of 
coordinate R and shows that, despite the similarities described above, the system is 
not a true quantum-mechanical simple harmonic oscillator. 

The uncertainty product ARAP follows from Eq. (1371) and is given by 



ARAP>±(±) . (38) 
At a fixed time, the minimum-uncertainty wave functions therefore satisfy 

{R-(R))1> = iA{P-(P)W , (39) 

where 

„ 2f(AR) 2 AR 

is a real constant. Defining the parameters 

/ / 1 

equation (]39f) becomes 

(9r-" + (7o + *Ti))V'(r)=0 . (42) 
The solution to Eq. (f42"D is 

^(r) = Nr a e- lor e-^ ir , (43) 

where iV is a normalization constant. These states form a three-parameter family of 
squeezed states minimizing the uncertainty relation ( |3"8"D at a fixed time. These are 
our radial squeezed states (RSS). 

Our derivation has demonstrated that, in any atom for which SQDT provides 
a good description, the RSS are candidates for Rydberg wave packets. Moreover, 
they have a relatively simple analytical form. In particular, these results apply to 
the alkali-metal atoms rubidium and potassium, for which the experiments on radial 
wave packets have been done to date. 

For the special choice a = I, corresponding to RSS with uncertainty ratio A equal 
to that of the ground state, the RSS specialize from the three-parameter family of 
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squeezed states to a two-parameter family of coherent states [|IJ. They then provide 
SQDT versions of the 'minimum-uncertainty coherent states' for three-dimensional 
systems that were introduced in Ref. flI5| . However, they correspond to states of 
Rydberg atoms with high angular momentum. This means they are unsatisfactory as 
descriptions of Rydberg atoms prepared by excitation with a single short laser pulse 
with no external field present, since these have large n and 1 = 1. We therefore do 
not discuss these coherent states further in the present work. 

D. Properties of Radial Squeezed States 

In this section, we present some time-independent properties of the RSS. The 
time-dependent behavior is discussed in Sec. V. 
First, we normalize ip{r). Imposing 

POO 

/ r 2 \ip{r)\ 2 dr = 1 (44) 
Jo 

gives 

N=^l . (45) 

This normalization requires a > — 1 and 70 > 0. 

Certain expectation values are of particular use in the description of the RSS and 
the comparison with other results. Here are a few key expectations, all calculated 
with normalized RSS. We find 

= 2^+3 ( 1 } = J» (46) 
270 r a + l 

(q + 2)(2a + 3) 1 2 7 q 2 

{ ' 2 7o 2 ' V ; (a + l)(2« + l) ' 1 ' 

(Pr> = -7i , (^ 2 ) = ^XT+7i 2 • (48) 

2a + 1 

It follows that the expectation value of the hamiltonian, which is the energy expec- 
tation value of the RSS, is given by 

< g> - 2( a + ifc + 1) N' ,(r +1)+a+1 )-^h +7 i ■ < 49) 
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From these expressions, it is apparent that the parameter a is further constrained 
by the requirement that the kinetic and potential energies of the RSS be separately 
normalizable, so that a> —\. 

It is also of interest to determine the uncertainties in the conventional coordinates 
r and p r . We find 

Al . = ^fH , = . (50) 

27o V2a + 1 

This implies that the RSS obey the uncertainty relation 



***=*v!tt ■ (51) 

As might be expected, the RSS are not minimum-uncertainty states in the usual 
variables r and p r , although by construction they satisfy the equality in the minimum- 
uncertainty relation (j38[). 

The possible value of the uncertainty product ranges from an arbitrarily large value 
when a is near its minimum of — |, through ArAp r = V3/2 ~ 0.866 when a = 0, 
to values falling as ArAp r ~ |(1 + l/a) for large a. As is discussed in Sec. V, large 
values of a are required to match key features of the RSS with the experimentally 
produced Rydberg wave packets at the first pass through the outer apsidal point of 
the orbit. We therefore expect the appropriate RSS to have uncertainty product close 
to I at that time. This feature agrees with the behavior of Rydberg atoms excited by 
a short laser pulse. Numerical calculations performed for hydrogen show that ArAp r 
is close to its minimum when the electron is near the outer turning point [III [|. 

The radial probability distribution, 

f(r) = r 2 \^(r)\ 2 = N 2 r 2a+2 e~ 2 ^ r , (52) 

determines the shape of the RSS. The function /(r) is asymmetrical in r. Its central 
maximum is located at r = ro, where 

a + 1 1 . , 

7o (-) 

The curvature of the envelope at r can conveniently be measured by the ratio 

C^tM.-lf . (54 , 
f(r ) a + 1 
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The specification of the RSS distribution is complete if Tq and Cq are known. 

The RSS envelope f(r) is a gamma distribution. The mean is just (r), while the 



the asymmetry of the function f(r), such as the momental skewness or the kurtosis, 
can be obtained as simple expressions in terms of a and 70. The detailed shape of 
the distribution is determined by its moments about the mean, defined as 



These can be found by direct integration or as the coefficients of the power-series ex- 
pansion of a relatively simple moment-generating function. In principle, the moments 
rrik provide a basis for a detailed comparison of the RSS with wave packets obtained 
from other theory or experiment. We do not further pursue this line of inquiry here. 

V. TIME EVOLUTION OF RADIAL SQUEEZED STATES 

At this stage, we have identified the RSS as candidate Rydberg wave packets and 
have discussed several time-independent features. In this section, we discuss some 
time-dependent properties, in particular the RSS time evolution and autocorrelation 
functions. We also compare our results for the RSS with results presented elsewhere 
in the literature. 

A. Time Evolution: Theory 

Since the RSS have a relatively simple form, their behavior as a function of time 
can be studied both numerically and analytically. We have done both. Our numer- 
ical integrations are performed using the Crank-Nicholson method. Our analytical 
studies apply to the time evolution of RSS and to their autocorrelation functions, 
including issues such as characterizing full and fractional revivals as a function of the 
supersymmetry integer 1(1) and the quantum defect 5(1). This subsection presents 
some analysis for time-dependent properties. In what follows, we define the origin 
t = to be the time that the packet first reaches the outer apsidal point. 

The full three-dimensional time-evolved wave packet can be written as a sum over 



variance is (Ar) 



2 . These are given in Eqs. ( ^6|) and ([)(]). Other standard measures of 




(55) 
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the complete set of SQDT states: 



y(r,t)=Y 10 (9,<P)Y,c n R n *i*(r)e- lE ^ , (56) 

n 

where I* — 1 — 5(1) + 1(1) is specified for p states. In principle, the expansion 
includes the continuum states. However, for cases of interest the contributions from 
continuum states are negligible. The coefficients c n are determined by requiring that 
the initial wave function \l/(r, 0) has radial piece given by the RSS ip(r) of Eq. fl43|) , 
i.e., \l/(r, 0) = Y w (9, <p)ip(r). This implies 

c n = (Rn*i*(r) I ^(r)> . (57) 

Using the expression for R n *i*(r) in Eq. (|j) and performing the integration, we get 
after some calculation the result 



N 



i 

Y(a + 1* + 3) (Y(n* + I* + 1) \ 2 n <^) 



T(l* + 1) V r(n*-Z*) ; [ n *( 7o + , 7l ) + i]^+ 1 



x 2 Fi (r + 1 - n *,a + r + 3;2r + 2;— 2 -— — -J . (58) 

V n (lo + *7i) + 1/ 

Here, 2-^1 is a hypergeometric function, and is the normalization constant given in 
Eq. fH). 

The coefficients in Eq. (^) can also be found via a suitable integration over the 
analytical propagator for the Coulomb problem obtained recently in Ref . |1T[ . For the 
examples of interest, we find that the coefficients in the sum over states in Eq. (|56| ) are 
strongly peaked around a central value n of n. Accurate numerical approximations 
can be obtained by truncating the sum to a subset of states centered on n. 

Another quantity of interest for the characterization of the time evolution is the 
absolute square of the autocorrelation function, A(t) = \ (*&(t) \ ^(0)) | 2 . In terms of 
the coefficients in Eq. (|5%D, it is given as 



A(t) = I E I c « I e~ lEntt ■ (59) 

n 

Studies of the behavior of a hydrogenic radial packet as a function of time show 
that wavefunction dispersion is followed by the appearance of revivals [1], |TTJ, |P2j| . We 



18 



have repeated these analyses for the RSS within the context of the SQDT. As the RSS 
evolves in time from its initial configuration at t = 0, it disperses. However, it is a 
trapped object, being bound to the atomic core. The dispersion therefore eventually 
leads to a situation in which what were initially the front and back parts of the packet 
interfere. This provides one signal for the eventual decoherence of the packet. For 
an RSS with eigenstate decomposition peaked about a principal quantum number n* 
and dominantly contained within a range 5n*, we find that in the SQDT that the 
interference time t* nt at which the packet effectively collapses is given approximately 
by the expression 

n* 

C = fe 23 ' (60) 

where T c * is given by Eq. (f24|). 

At a time substantially later than t- m t, the RSS reforms approximately into its 
original shape. For the SQDT, we find that this behavior occurs near a revival time 
t* iev given by 

C = K T d ■ (61) 

The new wave packet oscillates with periodicity equal to the classical orbital period 
T c *. We refer to this packet SIS db full revival 114 31. 

At certain times between t* nt and t* ev , the RSS gathers into r spatially separated 
packets called fractional revivals. A full analysis in the SQDT, along the lines of that 



presented in Ref. [11], suggests that these times are given by rational fractions 2p/q 
of the revival time t* ev , where p and q are relatively prime and where the number of 
packets is r = q if q is odd or r = q/2 if q is even. For the purposes of the discussions 
in the subsections below, we single out the particular fractional-revival times 

K = Kv ■ (62) 
The wave-function period of these fractional revivals in the SQDT is given by 

T; = ±T C * . (63) 
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With our conventions, the full revival corresponds to the case r — 1. 



B. Initialization 

To compare the properties of the RSS for different alkali-metal atoms with other 
theory and experiment, we need a procedure for determining the values of the pa- 
rameters a, 70, and 71 in a given situation. In a Rydberg atom that has been excited 
by a short laser pulse, the uncertainty product ArAp r is expected to be lowest at 
the first pass through the outer turning point |[], We therefore take this as our 
initialization point, and, as above, define it as the time origin t = 0. 

Fixing the initial form of the RSS requires the specification of three quantities. 
We choose these as the physical quantities (p r ), (r), and (H). The natural choice for 
the expectation value of the radial momentum at the apsidal point is zero. Similarly, 
a natural choice for the expectation value of r at the initialization point is the outer 
apsidal point of the orbit, 



*2 

n 



1 + 1 1 



+ 1) 



n* 2 



(64) 



Note that this differs from the outer apsidal point of the classical orbit, given in Eq. 
(Effi), by the use of the quantum-mechanical eigenvalue for the angular-momentum 
operator. Finally, one natural choice for the energy expectation value is the energy 
E n * = —l/2n* 2 of the central state n in the range of states excited by the short laser 
pulse. 

Thus, we initialize the RSS by choosing the parameters a, 70, 71 to satisfy the 
following equations: 

(Pr)=0 , (65) 

(r) = r * ut , (66) 
(H) = E n . . (67) 



The expectation values appearing in these equations are given in Eqs. (46), (48), and 



|). The first condition, Eq. (|65l) , can be inverted immediately to give 71 = 0. This 
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leaves two conditions to fix the two remaining parameters a and 70 in terms of the 
two quantities n* and I*. These in turn depend on the two quantities 1(1) and 8(1). 

For the Rydberg wave packets obtained in experiment the value of I* is much 
smaller than the value of n*. This means that the value of a is much greater than 
that of 7q. In the cases we consider below, a is large and exceeds 7q by more than 
six orders of magnitude. From Eqs. ( |5UD it follows that the initial squeezing in p r 
exceeds that in r by about six orders of magnitude: Ar/Ap r ~ a/7o ~ 10 6 . For 
more details about the uncertainty properties of the RSS, including a discussion of 



the time-dependence of the uncertainty product and ratio, see Ref. [^4 

C. Examples: Time Evolution 

In this subsection, we discuss two examples of the time evolution of RSS. The first 
is an RSS for hydrogen, with dominant energy component corresponding to principal 
quantum number n = 85. This example serves as a basis for comparison with other 
theoretical results. The second example is for rubidium. To determine the effects of 
the quantum defect on the time evolution, we choose for this case the same value of 
n, which corresponds to n* ~ 82.35. 

For our first example, the SQDT reduces to the usual Coulomb case for p states: 



1(1) = 5(1) = 0, /* = I = 1, and n* = n = 85. From the conditions ( |66"D and (67]), we 
determine that the RSS parameters are a ~ 168.225 and 70 — 0.0117465. This gives 
an initial uncertainty product of ArAp r ~ 0.50148. We find that initially the wave 
packet moves toward the inner turning point. The packet decreases in amplitude, 
broadens, and starts to oscillate as it approaches the origin. After the bounce, it 
regathers into a coherent packet as it approaches the outer turning point again. We 
find that the periodic motion of the RSS is consistent with the classical orbital period 
T* = T cl = 2vm 3 ~ 93.3 psec. 

Figures 3 and 4 show the envelope of the RSS at different times during the first 
orbit or so. In Fig. 3, the envelope of the RSS at t = is centered about the outer 
turning point near r out ~ 2n 2 ~ 14450 a.u. The direction of propagation is towards 
the origin. The envelopes at times t ~ 19.4 and 38.7 psec are also shown. Fig. 4 
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shows the RSS envelopes after reflection off the core, while approaching and passing 
through the outer turning point again. From left to right, the graphs are at times 
t ~ 58.1, 77.4 and 96.8 psec. The latter plot displays a decreased amplitude because 
the packet is a few picoseconds into its second orbit and is therefore moving towards 
the origin again. 

The generic behavior we find is consistent with the results of our perturbative 
calculation shown in Fig. 1. However, the perturbative calculation starts to break 
down after about one orbit, and comparisons cannot be made for long times. Our 
results compare well with Fig. 1(b) of Ref. [Q, which shows the Rydberg wave packet 
during the second part of its orbit, obtained by numerically integrating the time- 
dependent Schrodinger equation for an electron in hydrogen subject to a short laser 
pulse. This comparison suggests that the RSS are good models for radial wave packets 
in hydrogen. 

Our second example involves rubidium p states with n = 85, which have 5(1) ~ 
2.65, 1(1) — 3, I* ~ 1.35, and n* ~ 82.35. These values determine the RSS parameters 
as a ~ 162.91 and 70 — 0.0121233, with an initial uncertainty product ArAp r ~ 
0.50153. 

Figures 5 and 6 show the time evolution of the envelope of this RSS for rubidium. 
The time intervals in Figs. 5 and 6 are the same as those used for hydrogen in Figs. 3 
and 4. Some differences between the rubidium and hydrogen examples are apparent. 
These may largely be attributed to the quantum defect. For example, the outer 
apsidal point for rubidium is r out ~ 2n* 2 ~ 13560 a.u. and is closer to the core than 
the corresponding point for hydrogen. Also, comparison of Figs. 3 and 5 shows that 
at t ~ 38.7 psec the rubidium RSS is oscillating significantly more than the hydrogen 
one. This is a consequence of the reduction of the classical orbital period of rubidium, 
T c * = 27m* 3 ~ 84.9 psec, relative to that of hydrogen. At 38.7 psec, the rubidium RSS 
is closer to halfway through its orbit and hence exhibits more quantum-like behavior. 

D. Examples: Autocorrelation Functions 

In this subsection, we provide some examples of autocorrelation functions for 
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RSS. Three cases are considered. The first two use the same RSS for hydrogen and 
rubidium considered in the previous subsection. The third example we use is an RSS 
for potassium with n = 67, which is convenient for comparison with experiment. 

Our first example is the RSS for hydrogen with n = 85. The associated parameters 
a and 70 are given in the previous subsection. The discussion in Sec. V A suggests 
that the qualitative features of the motion are governed by the classical orbital period 
T c i ~ 93.3 psec, by the interference time t int ~ 4T c i obtained with a spread 5n ~ 7, 
and by the revival time t rev ~ 2.6 nsec. 

Figure 7 shows the absolute square of the autocorrelation function for hydrogen 
with n = 85. Near t = the autocorrelation exhibits peaks spaced by the classical 
orbital period, reflecting the oscillatory motion between the outer and inner apsidal 
points expected at early times. About four orbits can be distinguished before the 
packet collapses. At times near the predicted value of t Tev , peaks spaced by the 
classical orbital period reappear. The figure also reveals the appearance of fractional 
revivals corresponding to the values r = 4, 3, and 2 at times ~ 0.65 nsec, £3 ~ 0.88 
nsec, and ti ~ 1.3 nsec, in agreement with the theory. As expected, these fractional 
revivals exhibit wave-function periodicities of T 4 ~ 23.3 psec, T 3 ~ 31.1 psec, and 
T 2 ~ 46.7 psec, respectively. Note that the maximum values of A near the times 
t = t r decrease with increasing r, as might be expected from the spatial separations 
and numbers of the fractional revivals. 

Our next example is the RSS for rubidium with n = 85. As described in the 
previous subsection, this corresponds to n* ~ 82.35. The SQDT predicts a modified 
behavior of this RSS, with key features determined by the classical orbital period 
T c * ~ 84.9 psec, by the interference time t* nt ~ 4T C *, and by the revival time t* ev ~ 2.3 
nsec. Note that although the interference time is different from that in hydrogen, it 
remains given by about four classical orbital periods because the values of n and Sn 
do not change significantly; cf. Eq. (j60|) . 

Figure 8 shows the absolute square of the autocorrelation function for this exam- 
ple. The overall behavior is similar to that of hydrogen. In agreement with the pre- 
dicted interference time, about four orbits occur before wave-packet collapse. Larger 
correlations separated by the classical orbital period reappear near the theoretical 
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value of t* ev . Fractional revivals with r = 4, 3, and 2 can be seen near t\ ~ 0.58 
nsec, t| ~ 0.78 nsec, and t\ ~ 1.2 nsec, as expected. The associated periodicities are 
T 4 * ~ 21.2 psec, T 3 * ~ 28.3 psec, and T* ~ 42.5 psec. 

Our last example is for the RSS in potassium with n = 67. The associated p 
states have 5(1) ~ 1.71, 1(1) = 2, I* ~ 1.29, and n* ~ 65.29. The corresponding RSS 
parameters are a ~ 120.306 and 70 — 0.0142822. The initial uncertainty product is 
ArAp r ~ 0.50207. For this case, the classical orbital period is T* x ~ 42.3 psec. With 
5n* ~ 7 still, the interference time is reduced to t* nt ~ 3T C *, i.e., about three orbits are 
expected to occur before the packet collapses. The revival time is now significantly 
smaller, t* ev ~ 0.92 nsec. 

The absolute square of the autocorrelation function for this potassium RSS is 
shown in Fig. 9. Once again, the features agree with the theory. Only the fractional 
revivals with r = 3 and 2 can be clearly seen, however. They occur near t$ ~ 0.31 
nsec and t\ ~ 0.46 nsec, as predicted. The wave-function periodicities are T 3 * ~ 14.1 
psec and T 2 * ~ 21.2 psec, respectively. 

This example can be compared to the experiment presented in Ref. 0, in which 
potassium radial wave packets were excited in a pump-probe experiment. Oscillations 
in the photoion signal were observed that are consistent with T c * and t* ev calculated 
from SQDT for a potassium RSS with n* ~ 65. The revival peaks observed would 
correspond to the autocorrelation peaks exhibited in Fig. 9 starting at about 0.8 nsec. 
Similarly, examination of the experimental results presented in Refs. |], [| concerning 
the observation of fractional revivals with r = 2 reveals behavior consistent with the 
analysis above. 

VI. SUMMARY AND CONCLUSIONS 

In this paper, we have investigated two possibilites for describing Rydberg wave 
packets obtained by atomic excitation with short laser pulses. The approaches are 
somewhat complementary and both are developed in the context of SQDT, which 
can be viewed as an analytical model for single-body effective wave functions of 
Rydberg atoms. As such, the results we obtain are relevant both for hydrogen and 
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for non-hydrogenic systems, including in particular the alkali-metal atoms used in 
experiments. 

The first approach is a direct extension of time-dependent perturbation theory 
for the generation and evolution of the wave packets. The success of this approach 
suggests that the SQDT may be of use in calculations of this type carried out for 
alkali-metal atoms. 

We have concentrated primarily on the second approach, which is an attempt to 
describe Rydberg wave packets in terms of a class of analytical squeezed states, the 
RSS. Using the SQDT, we succeeded in constructing the RSS as a three-parameter 
family of analytical states minimizing an uncertainty relation in nonstandard coordi- 
nates. The analytical form of the RSS, which are gamma distributions, means that 
relatively simple and/or analytical expressions can be obtained for various proper- 
ties. The use of SQDT ensures that the RSS provide candidates for the analytical 
description of Rydberg wave packets in atoms other than hydrogen. 

In terms of the radial coordinates r and p r , the RSS have close to minimum un- 
certainty for the cases we consider. To match with other theory and experiment in 
specific cases, we initialize the RSS at the outer apsidal point of the orbit. The three 
parameters are fixed by matching the expectation values of the position, momentum, 
and energy at that point. The ensuing behavior is completely determined. We derive 
relatively simple expressions for the time-evolved RSS and for their autocorrelation 
functions. We also predict the dependence of characteristic features of the behav- 
ior on the quantum defects. These include the initial orbital period, the time of 
packet collapse, and the full and fractional revival times and their associated orbital 
periodicities. 

We have presented examples of time evolutions and autocorrelation functions both 
for hydrogen and for alkali-metal atoms. The behavior agrees with our theoretical 
expectations. The properties and behavior of the RSS are also in agreement with 
results on radial probability distributions, time evolution, and autocorrelation func- 
tions obtained by other authors for Rydberg wave packets. It appears that the RSS 
provide an analytical tool for describing and distinguishing experimental and theo- 
retical features of Rydberg-atom radial wave packets produced by short-pulsed laser 
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fields. 
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Figure Captions 

Fig. 1: Hydrogen radial wave packet with n = 85. The unnormalized radial proba- 
bility distribution f(r) is plotted as a function of radial distance r in a.u. From 
left to right, the graphs are shown for the times t = |T c i, |T cl , |T cl , and |T cl , 
where T cl = 93.3 psec. 

Fig. 2: Rubidium radial wave packet with n = 85. The unnormalized radial proba- 
bility distribution f(r) is plotted as a function of radial distance r in a.u. From 
left to right, the graphs are shown at the same times as in Fig. 1: t — 10.4, 
20.7, 31.1, and 41.5 psec. 

Fig. 3: Hydrogen RSS with n = 85: first part of the orbital motion. The unnor- 
malized radial probability distribution f(r) is plotted as a function of radial 
distance r in a.u. From right to left, the graphs are shown at the times t — 0, 
19.4, and 38.7 psec. 

Fig. 4: Hydrogen RSS with n = 85: second part of the orbital motion. The un- 
normalized radial probability distribution f(r) is plotted as a function of radial 
distance r in a.u. From left to right, the graphs are shown at the times t = 58.1, 
77.4, and 96.8 psec. 

Fig. 5: Rubidium RSS with n = 85: first part of the orbital motion. The unnor- 
malized radial probability distribution /(r) is plotted as a function of radial 
distance r in a.u. From right to left, the graphs are shown at the same times as 
in Fig. 3. 

Fig. 6: Rubidium RSS with n = 85: second part of the orbital motion. The un- 
normalized radial probability distribution /(r) is plotted as a function of radial 
distance r in a.u. From left to right, the graphs are shown at the same times as 
in Fig. 4. 

Fig. 7: Absolute square of the autocorrelation function for hydrogen RSS with n = 
85 as a function of time in nanoseconds. 

Fig. 8: Absolute square of the autocorrelation function for rubidium RSS with n = 
85 as a function of time in nanoseconds. 

Fig. 9: Absolute square of the autocorrelation function for potassium RSS with n = 
67 as a function of time in nanoseconds. 
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